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Abstract 

The unquenched spectral density of the Dirac operator at / is complex 
- - - , and has oscillations with a period inversely proportional to the volume and an 

^ I amplitude that grows exponentially with the volume. Here we show how the 

O ' oscillations lead to the discontinuity of the chiral condensate. 

<N : 

. 1 Lessgo ! 

O ■ The sign problem in QCD at non-zero baryon chemical potential, p,, has been an obstacle for more than 

OO . two decades. New innovations have given exiting results (D even though the core of the sign problem in 

^ \ QCD remains. Here we address a closely related and equally long-standing problem (see e.g. Il2l): 

^ ' How does chtal symmetry breaking at /i / manifest itself in the spectrum of the Dirac operator? 

00 : 

T— I ■ At /i = the Dirac spectrum is located on the imaginary axis and the discontinuity of the chiral con- 

^ ! densate at zero quark mass, m, is proportional to the eigenvalue density at the origin (the Banks Casher 

^ I relation). While this is quite intuitive the situation for /u 7^ is quite puzzling. The chemical potential 

■ sends the eigenvalues, z^, of D + off into the complex plane while the discontinuity of the chiral 
^ . condensate remains. The solution 13| of this (Silver Blaze [4J) puzzle is in a sense the sign problem: 

\ Due to the sign problem the eigenvalue density becomes a complex and strongly oscillating function and 

■ the oscillations lead to the discontinuity of the chiral condensate. To show this we start from the exact 
solution |5 1 of the eigenvalue density for ^ ^ Aqcd and ^ l/V^ ((...) is the quenched average) 

I . (Ek + iy- Zk) det^/(D + + m)) 

. (det f{D + ijjo + m)) 
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and compute the chiral condensate S = limm^o.v^oo '^Nfi'm) using 

Sat (m) = — / dxdy — — ■ . (2) 

Given the nature of the sign problem it is perhaps not surprising that the oscillations of the eigenvalue 
density Q have a period inversely proportional to the volume and an amplitude growing exponentially 
large with the volume The results El|5||^ for the eigenvalue density are ideally suited to resolve these 
oscillations since they describe eigenvalues \z\ ~ 0{l/{TiV)). In [5'| the eigenvalue density was derived 
using non hermitian random matrix theory 1 8 1 while in [61 it was derived from the chiral Lagrangian 
using the replica method |9 1. 

On the next pages we explain the original |3]| direct computation of S n^. (m) through several figures. 
2 The phase quenched way 

Before going to the real problem we take a short aside and look at the way chiral symmetry breaking 
affects the spectral density of the Dirac operator in phase quenched QCD (tantamount to QCD at non- 
zero isospin chemical potential). This case serves to show how the perhaps more familiar mean field 
results are a special limit of the exact results. The real and positive eigenvalue density 

(Efc + W- ^k) |det(Z) + /i7o + 

P2n{x, y, m; p) = , , \ y;2^y (3) 

(|det(Z; + /U7o + m)|"=") 
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Fig. 1: Top: The eigenvalue density in tiie complex eigenvalue plane with two phase quenched flavors for fiF\/V = 8 and 
mEV = 80. Bottom: The corresponding chiral condensate as a function of mT,V, again for fiFVV = 8. 

is plotted in the top panel of figure ^ (the explicit expression is given in (76) of |6|). As expected 
the eigenvalues have spread out away from the imaginary axis. Notice the absence of eigenvalues at 
z = ibm and at z = 0. Using © the corresponding chiral condensate follows, see the lower panel of 
figure n The drop of the chiral condensate as the quark mass enters the eigenvalue distribution is easy 
to understand by an electrostatic analogy and is consistent with lattice measurements flQl. Furthermore, 
taking mSF > 1 and /i^fV > 1 the mean field result = mT? / {2^? F"^) for m^/(2/i) < 1 and 

S otherwise) is reproduced by the exact result. {F is the pion decay constant.) 

3 The unquenched way 

In full QCD the eigenvalues of the Dirac operator for a given gauge configuration are indistinguishable 
from what we would have obtained in the phase quenched case. However, due to the phase of the 
fermion determinant, the average spectral density is entirely different. To show this we have plotted the 
unquenched eigenvalue density in the top panel of figure|2l Its explicit expression is given in (73) of [6|. 
Here it is sufficient to note that this expression naturally separates into the quenched spectral density and 
a remainder 

PNf{x,y,m;n) = pQ{x,y;fi) + pu{x,y,m; ^i). (4) 

These two parts are shown in the two lower panels of figure |2l The quenched eigenvalue density is 
real and positive and behaves as the phase quenched eigenvalue density (except for the dip at the quark 
mass). The unquenched part contains the complex oscillations of which we have only displayed the real 
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Fig. 2: The eigenvalue density of the full QCD Dirac operator for one flavor (real part only) and fiFVV = 8 and mT,V = 80. 
The bottom figure shows the difference between the full (top) and quenched (middle) eigenvalue density. Notice the similarity 
between quenched and the phase quenched spectral density shown in figureQ For y = the oscillations start |a;| = m. 

part. For the values of the parameters mSF and fi'^F'^V given in figure |2j the maximum value of the 
amplitude of the oscillations is 400 times larger than the plateau of the quenched spectral density. For 
better illustration the oscillation have been clipped. Inserting the density Q into ^ leads to two terms 



(5) 



1.2 r 

0.9 : 

0.6 : 

0.3 : 

: 

-0.3 : 

-0.6 : 

-0.9 _ 
-1.2 



-160 -120 



-80 



-40 



40 



80 120 160 




-40 

mLV 



80 120 160 



Fig. 3: The chiral condensate as a function of quark mass obtained from the three densities in figure|2|using J2}. 
Top: The chiral condensate in full QCD with the discontinuity at m = (independent of /i) is the sum of the two below. 
Middle: The quenched chiral condensate (dependent on /i) 1 13 1. No discontinuity at m = 0. 
Bottom: The contribution from the oscillating part of the eigenvalue density (dependent on /i). 



The quenched part, Eg, drops to zero as m comes inside the support the eigenvalue density. However, 
the unquenched oscillating part exactly makes up for this and leaves the full chiral condensate S jv^ (m) 
independent of fi and therefore equal to the result for /i = 1 11 1 (see figure |3ll. While Sq is built up 
by the eigenvalue density inside the quark mass (the contributions from |x| > m cancel each other), the 
unquenched part S(/ is built up by the density pu outside the quark mass (not easy to understand by an 
electrostatic analogy). This is possible since the oscillations have a period of order 1 /V and an amplitude 
that grows exponentially with V LIJ. Moreover, the contribution from the unquenched part comes from 
the entire oscillating region and not from the boundary of the support of pu. 
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4 Conclusions 

While the discontinuity of the chiral condensate is usually associated with a dense spectrum of Dirac 
eigenvalues on the imaginary axis near zero, the sign problem allows for an alternative mechanism: The 
discontinuity of the chiral condensate arises from strong oscillations of the eigenvalue density. Having 
a period of order 1/V and an amplitude that grows exponentially large with V the oscillations directly 
reflect the sign problem. The oscillations are present as soon as the quark mass hits the support of the 
eigenvalue density, i.e. when the phase quenched theory enters the BEC phase. A lattice simulation in 
this region therefore must be able to deal with the oscillations in order to study the correct mechanism 
for spontaneous chiral symmetry breaking in QCD Ill2i . 
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